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This paper reports a study of B meson oscillations using hadronic Z decays with 
two identified leptons, and updates a previous publication by including data collected in 
1994. Decay times are reconstructed for each of the semileptonic B decays by forming 
vertices which include the lepton and by estimating the B meson momentum. The mass 
difference, Amj, between the two mass eigenstates in the B^ system is measured to 
be 0.430 ± 0.043 _o'o30 P s > where the first error is statistical and the second error is 
systematic. For the Bg system, a lower limit of Am s > 2.2 ps^ 1 is obtained at 95% C.L. 
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1 Introduction 



In the Standard Model, a second-order weak transition transforms neutral B mesons into their 
antiparticles JO. The neutral B mesons therefore oscillate between particle and antiparticle 
states before decaying. The frequency of the oscillation depends on the top quark mass, the 
Cabibbo-Kobayashi-Maskawa matrix elements, and meson decay constants. By analogy with 
the K° case and neglecting CP violation, the mass eigenstates, |Bi) and |B 2 ), of B° (q=d or s) 
can be described as follows: 

|Bi) 
|B 2 ) 

If a B° is produced at time t = 0, the probabilities of having a B° or a B° at proper time t areQ 

P B o(t) = -e-* cos 2 (-J^) 
Pgo(t) = -e-^ sm 2 

where r is the B° lifetime. The frequency of the oscillation is given by Am q , the mass differ- 
ence of the two mass eigenstates (Am q = — m^ 2 ). For B° — B^ mixing, time-integrated 
measurements from ARGUS and CLEO give x& = Am d r = 0.67 ± 0.08 0, |3|. Published 
measurements of the frequency of B d — B d oscillations made at LEP are available using several 
different techniques [|], |5], H 1] • Lower limits on Am s have been reported by the ALEPH || j(| 
and the OPAL || |7| collaborations. 

Extracting information on CKM matrix elements from the measurements of Amd and Am s 
is prone to large uncertainties due to poorly known meson decay constants. These uncertain- 
ties can be reduced by considering the ratio Am s / Arrid- Given the present knowledge of V ts 
and Vtd one expects Am s to be of the order of 10 ps _1 [0. Using dilepton events in data 
collected between 1991 and 1993 ||, we studied B[J — B° and B° — B° oscillations, reporting 
Am d = 0.462 ±° Q ° f 3 t^l ps" 1 and Am s > 2.2 ps" 1 at 95% C.L. We update these results by 
including data collected in 1994. The technique is the same as that reported previously 0. 
Hadronic Z° decays with two lepton candidates, one in each thrust hemisphere, are selected. 
The reconstruction of a secondary vertex that includes the lepton is attempted for each lepton 
candidate, yielding an estimate of the decay length of the b hadron. This is combined with an 
estimate of the relativistic boost of the b hadron to give the proper decay time. The likelihood 
of each event is calculated as a function of Amj and Am s according to the measured proper 
times and the charge combination of the two leptons. The result for Amd and the lower limit 
on Am s are then obtained using a maximum likelihood technique. 



2 Event Selection and Simulation 



1 The contribution of Ar, the difference between the total decay widths of the mass eigenstates, to the 
oscillations is expected to be negligible and has been ignored. 



^(|b;> + ib;», 



-^(|b;> - ib;». 
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2.1 Event Selection 



The analysis is performed on data collected by OPAL in the vicinity of the Z° peak from 1991 



to 1994. The OPAL detector has been described elsewhere [IT, T2\. Hadronic Z° decays are 
selected using criteria described in [[13]. A cone jet algorithm [TJ]] is used to classify tracks and 
electromagnetic clusters not associated to tracks into jets. The size of the cone is chosen so 
as to include nearly all the decay products of a b hadron into one jet. The jets also include 
particles produced in the fragmentation process, which originate from the e + e~ collision point. 
A total of 2 874 660 hadronic events satisfy the event selection criteria. 

Electrons are identified using an artificial neural network |J which is trained on a sample 
of simulated hadronic Z° decays. Electrons from photon conversions are rejected as in JT3| . 
Muons are identified as in 



16 



Lepton candidates are required to satisfy p > 2.0 GeV and 
| cos#| < 0.9. Additional kinematic criteria are imposed to reduce the fraction of leptons in the 
sample coming from cascade decays of the type b — > c — > £. 

The techniques for secondary vertex reconstruction and proper time estimation are described 
in Pj. Dilepton events with at least one reconstructed vertex are selected. 



2.2 Event Simulation 



Monte Carlo events are used to predict the relative abundances and decay time distributions for 



lepton candidates from various physics processes. The JETSET 7.4 Monte Carlo program [17 

is used to generate Z° — > qq events which are 
The fragmentation of b and c quarks is 



19 



with parameters tuned to OPAL data 
processed by the detector simulation program 
parametrised using the fragmentation function of Peterson et al. 
hadrons given by the central values in Table |1|. 



20] , with (xe) for b and c 



Quantity 


Value 


(x E )b 


0.697 ±0.013 fl6| 


(xe)c 


0.51 ±0.02 I IS | 


B(b -> £) 


(10.5 ±0.6 ±0.5)% 16] 


B(b^c^ t) 


(7.7 ±0.4 ±0.7)% pfl 


B(b^c^ i) 


(1.3 ±0.5)% @ 


M(B°) 


5.48 GeV 


M(A b ) 


5.62 GeV 


r B +/r B o 


1.03 ±0.06 § 


T B°J T B° 


1.03 ±0.08 U 


r Ab/ r B° 


0.73 ±0.06 § 


<Tb> 


1.55 ± 0.02 ps @| 



Table 1: The parameters used for the Monte Carlo simulation. 



Standard Model values of the partial widths of the Z° into qq are used jpT] . The mixture of 
c-flavoured hadrons produced both in Z° — > cc events and in b hadron decays is as prescribed 
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Ill 



15]. The semileptonic branching ratios of charm hadrons and associated uncertainties are 



also those of |L5j. The central values in Table [l] are taken for the inclusive branching ratios 
for b — y £, b — > c — > £ and b — > c — > i. The semileptonic branching ratios of the individual b 
hadrons are assumed to be proportional to the lifetimes. The models used in describing the 
semileptonic decays of b and c hadrons are those used in determining the central values in 
15 1 . The asssumed masses for B° and Ab particles are also given in Table [l]. The lifetimes of 



b hadrons used in this analysis were taken from the world average values |j, as indicated in 
Table |. 



3 Fit Results for Am d 

The numbers of dilepton events with at least one secondary vertex constructed for the combi- 
nation of e-e, e-/i and are listed in Table ^], separately for like-sign and unlike-sign dilepton 
events. Also included is the total number of secondary vertices reconstructed in these events. 





e-e 


Grfl 




total 


total vertices 


unlike-sign 


891 


1791 


1070 


3752 


5971 


like-sign 


377 


780 


448 


1605 


2573 



Table 2: The numbers of dilepton events with at least one secondary vertex reconstructed 
for the combinations e-e, e-/x and separately for like-sign and unlike-sign leptons. Also 
indicated is the total number of secondary vertices reconstructed in unlike-sign and like-sign 
dilepton events. 



In order to study Amj and Am s , the likelihood of the event sample is calculated as a 
function of these parameters. The construction of the likelihood function follows the procedure 
described in the previous paper |J . The true proper time distribution is described by a physics 
function for each source of events. The B mixing is also described by the physics function. The 
reconstructed time distributions, /(£), are then obtained by convolving the physics function 
with resolution functions, P(t,t'), which describe the proper time resolution for each source. 
For example, in the absence of mixing, 

f(t)=- r e -TP(t,f)df 

r Jo 

for a source with lifetime r. The resolution functions are complicated functions, which must 
include a description of misreconstruction near t = 0, even when the true proper time is large. 
The probability of this misreconstruction depends on the true proper time, and the description 
of the resolution function was modified from the previous paper to describe this better. The 
resolution function has the following form: 

P{t, t') = C[(l - exp {--))v{t, + exp (--)u(*)] 

a a 

where t and t' are the reconstructed proper time and the true proper time, respectively. The 
functions v(t,t') and u{t) describe the reconstructed proper time distributions for the correctly 
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reconstructed and misreconstructed vertices respectively. The parameter C is a normalization 
factor, while a is a parameter to describe the dependence of the misreconstruction probability 
on the true proper time. Distributions of t and t—f are shown in Figure [I] for three slices of true 
proper time t' for Monte Carlo b decays. The fitted resolution function, which is superimposed 
in the figure, describes the distributions well. 

To determine Amd a three parameter fit is performed, varying Amj simultaneously with 
the cascade fraction, the fraction of lepton candidates in Z — > bb decays that are due to 
b — ► c — > £ decays, and the B° production fraction, the fraction of b quarks that give rise 
to B° mesons. Gaussian constraints reflecting the systematic errors on these two parameters 



are imposed. The relative uncertainty in the cascade fraction is taken to be ±15% fTS] , which 
includes uncertainties due to branching fractions, decay modelling and detector simulation. 
The B° production fraction, / s , is constrained both by direct measurements, giving a rate of 
(11.1 ±2.6)% || relative to all weakly decaying b hadrons, and by the measured average mixing 
rate of b hadrons, x = 0.126 ± 0.008 || together with knowledge of the equivalent parameters, 
Xd an d Xs, for B° and Bg mesons (xd(s) = 0.5 x x\, B J{l ± x\r s \)). This is equivalent to the 
constraint / s = (11.2lj;g)% || except that the values of Xd and Xs are calculated from the 
values of Am^ and Am s in the fit, together with the appropriate lifetimes. 

In the fit the Bg oscillation parameter is fixed at Am s = 10.0 ps" 1 . The result of the fit is 
Am^ = 0.430 ± 0.045 ps _1 . The fitted value of the cascade fraction is 0.078 ± 0.007 compared 
to the nominal value of 0.066, as calculated from the Monte Carlo sample. The fitted value of 
/ s is (12.7 ±1.9)%. 



Figure [2] shows the distribution of decay times for all leptons in the dilepton sample and 
separately for leptons in like-sign and unlike-sign events. The curves are the results of the 
likelihood fit. 

The fraction of like-sign leptons as a function of proper decay time, 

N LS (t) 



nit) 



N us (t) + N LS (t) 



is plotted in Figure || for data, where N LS (t) (N us (t)) is the reconstructed time distribution for 
leptons in like-sign (unlike-sign) events. In the figure, the expected curve for Amd = 0.430 ps _1 
is shown as the solid line. The fitted values of the cascade fraction and the fraction of leptons 
from Bg decays are used. Events in which vertices have been reconstructed in both thrust 
hemispheres enter the plot twice. 



4 Systematic Errors on Amd 

In the three parameter fit, the error on Amd is a combination of statistical error and systematic 
error due to the constraints on the cascade decay fraction and the Bg production fraction. The 
systematic error from the cascade decay fraction is estimated by repeating the fit with the 
central value of the constraint for the cascade fraction changed by ±15% or —15% (the sys- 
tematic uncertainty on this parameter) from its nominal value. The systematic error resulting 
from the Bg production fraction is obtained in a similar way. The statistical error on Amd is 
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±0.043 ps l , obtained by subtracting in quadrature these two systematic errors from the fit 
error. 

The uncertainty due to the resolution function description is assessed by repeating the 
parametrisation using Monte Carlo events in which the tracking resolution is degraded by 10% 
2"2" or improved by 10%. The uncertainty in the background from Z° — > cc events is taken 
to be ±30% due to uncertainties in the branching fractions and modelling of semileptonic 
charm decays, the relative production rates of charmed hadrons, and the uncertainty in the 
partial width for Z° — > cc. The production rates of B° and B + are assumed to be equal and 



the b-baryon production rate is assumed to lie in the range (9 ± 4)%. The fraction of D** 
produced in decays of b hadrons was assumed to be B(b — > D**) = 0.36 ± 0.10. Variations in 
the efficiency of the secondary vertex reconstruction as a function of decay length are found 
to have a negligible effect on Amj. Uncertainties in the source composition due to Monte 
Carlo statistics are negligible. The B lifetime variations are performed by changing the ratios 
of individual B lifetimes while keeping the average lifetime fixed at the LEP, CDF and SLD 
average value, 1.55 ps -1 . 

The summary of the sources and estimated values of systematic errors is given in Table |^. 
The sum of these systematic errors in quadrature is 5Am d =to;o3o P s 



Source of uncertainty and range 


5Ama ps 1 


cascade decay fraction (±15%) 


-0.010 ±0.011 


B° fraction (see text) 


-0.006 ±0.006 


resolution function (±10%) 


±0.010 -0.010 


lepton misidentification (e: ±30%, /i : ±20%) 


-0.001 ±0.000 


charm background (±30%) 


-0.006 ±0.001 


b-baryon fraction (±0.04) 


±0.010 -0.009 


B(b -> D**) (±0.10) 


±0.003 -0.003 


r B+ /r B o = 1.03 ±0.06 
rB o/r B o = 1.03 ±0.08 
ta1/t b o = 0.73 ± 0.06 


±0.019 -0.023 
-0.000 ±0.001 
±0.004 -0.004 


Am s = 2-20 ps" 1 


±0.004 -0.000 


Total systematic error 


±0.028 -0.030 



Table 3: Summary of systematic errors on the Am^ measurement. For each source of error, 
the first error quoted results from varying the parameter in the positive sense. 



5 Fit Results for Am s 

We use the dependence of the likelihood on the assumed value of Am s to derive a lower limit. 
To account for systematic errors when setting the limit, we produce a likelihood curve as a 
function of Am s that includes these effects. This is achieved by maximizing the likelihood with 
respect to the value of each relevant parameter, constrained by a Gaussian error corresponding 
to its uncertainty, at each value of Am s . The constraints are those shown in Table |3|. In 
addition, Am d is treated as a systematic uncertainty constrained by the average Ama from 
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analyses using reconstructed D* mesons: Am d = 0.52 ± 0.05 ps" 1 @, §. The exception to 
this scheme is the treatment of the resolution function description. In this case, three curves of 
In £ were calculated, including all other systematic uncertainties, assuming the default tracking 
resolution or assuming the tracking resolution was degraded or improved by 10%. The smallest 
of the three values of — Aln£ was taken at each value of Am s . The solid curve in Figure |] 
shows the difference in log-likelihood from the maximum as a function of Am s with systematic 
errors included. 

We set the limit on the basis of the difference in log-likelihood, A In C, with respect to the 
maximum value. A Monte Carlo technique was used to determine the correspondence between 
confidence levels and values of A In £ as a function of Am s . This approach is found to be more 
reliable than the approach of the previous paper ||, where — Aln£ = 1.92 was assumed to 
correspond to 95% confidence level. Many data sets, of the same size as the real data sample, 
were simulated using a fast Monte Carlo and fitted in a manner similar to the data. The 
main systematic errors affecting the Am s result were simulated by allowing the parameters of 
the Monte Carlo to vary independently for each data set. The corresponding parameters were 
allowed to vary under Gaussian constraints in the fit. The exception to this was the parameter 
governing the proper time resolution, which was treated in the same way as in the data fit. For 
each simulated data sample a single value of A ln£ = ln£ max — ln£(Am*) was extracted, where 
Am* is the generated value of Am s . The Aln£ corresponding to 95% confidence was defined 
to be that value above which lay only 5% of the simulated data samples. This procedure was 
performed for input values of Am* = 1.0,2.0,4.0 and 8.0 ps _1 , using 3000 Monte Carlo data 
sets at each value of Am*. The results of this study are shown as the dashed line in Figure 
We exclude the region of Am s < 2.2 ps _1 at 95% C.L. 

To assess the importance of the systematic errors, we studied the log-likelihood as a func- 
tion of Am s , while fixing all other parameters. These parameters were set to the values that 
maximized the log-likelihood at the preferred value of Am s in the procedure described above. 
The result is shown as the dotted curve in Figure |j. We conclude that the systematic errors 
have only a minor effect on our result. 

Using the data sets simulated with the fast Monte Carlo referred to above, we were able 
to check the analysis technique and study the expected sensitivity to Am s . The results of 
these studies are shown in Figure ||. Each row corresponds to a different generated value of 
Am s , indicated by the Am* at the right edge of the plot. The left column shows the fitted 
value of Am s for each trial. In the right column are normalised cumulative distributions of 
the log-likelihood difference between the value at the fitted maximum and the value at the 
generated Am s for each trial. The sensitivity of this analysis is good for Am* < 4 ps _1 , but is 
lost between 4 ps _1 and 8 ps^ 1 . 



6 Conclusion 

We have measured the oscillation frequency Am d by measuring the proper time of B meson 
decays and tagging the charges of leptons in both thrust hemispheres. The B° — B° oscillation 
parameter is measured to be: 

Am d = 0.430 ± 0.043 t° f ps" 1 , 
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corresponding to (2.83 ± 0.28±q;2o) x 10 4 eV. This result is consistent with and supersedes the 
result using 1991-1993 data. 

The Am d value is consistent with the OPAL results Am d = 0.548±0.050 ±0:019 P s_1 from 
data containing D* ± mesons and leptons, and Amj = 0.444 ±0.029 ±0 017 ps _1 [|7| from inclusive 
lepton events. Combining these results, taking into account correlations in the systematic errors, 
we find 

Am d = 0.467 ± 0.022 ps" 1 . 

The small statistical correlations between the results were found to have a negligible effect. 
This result is consistent with previous measurements J|, §]. Using r B o = 1.56 ± 0.06 ps, the 
combined OPAL value gives x d = 0.73±0.04±0.03, where the last error is due to the uncertainty 
in r B o. This value is also consistent with the average of ARGUS and CLEO measurements, 
x d = 0.67 ±0.08 g §. 

We obtain a lower limit on Am s at 95% confidence level: Am s > 2.2 ps^ 1 . This limit is less 
constraining than the ALEPH results||, [| and a recent OPAL result 0. 
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Figure 1: The distributions of reconstructed proper time, t, and t — t' in three slices of the 
true proper time t' for leptons from primary b hadron decays in the Monte Carlo. Also shown 
is the parametrisation of these distributions. 
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Figure 2: The proper time distributions for all leptons in dilepton events (top) for which a 
vertex is found, and for those leptons in like-sign (centre) and unlike-sign (bottom) events. The 
curves represent the results of the maximum likelihood fit. 
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Figure 3: The fraction of like-sign leptons as a function of proper decay time: lZ(t). The solid 
curve represents the expectation with Am d set to 0.430 ps _1 and Am s set to 10.0 ps _1 . 
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Figure 4: The difference in log-likelihood from the maximum value is shown as a function of 
Am s . The solid curve includes the effect of systematic errors, while the dotted curve includes 
only statistical errors. The dashed curve shows the 95% C.L. 
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Figure 5: The results of fits to 3000 toy Monte Carlo datasets are shown. The Am* value 
indicates the generated value of Am s in ps _1 . The left-hand column shows the distribution of 
fitted values of Am s for each Am* value. The right-hand column shows the normalised cumu- 
lative distribution of the difference between the log-likelihood values at the fitted maximum 
and at the generated value. The arrow indicates the 95% confidence level value. 
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